TO APPEAR IN COSMOS ApJS SPECIAL ISSUE 2007 
Preprint typeset using IAT^X style cmulatcapj v. 10/09/06 



THE XMM-NEWTON WIDE-FIELD SURVEY IN THE COSMOS FIELD: VI. STATISTICAL PROPERTIES OF 

CLUSTERS OF GALAXIES/ 

A. Finoguenov 1-2 , L. Guzzo 3 , G. Hasinger 1 , N. Z. Scoville 4 - 5 , H. Aussel 6 , H. Bohringer 1 , M. Brusa 1 , P. Capak 4 , 
N. Cappelluti 1 , A. Comastri 7 , S. Giodini 3 , R. E. Griffiths 13 , C. Impey 8 , A. M. Koekemoer 9 , J. -P. Kneib 10 , A. 
Leauthaud 10 , O. Le Fevre 10 , S. Lilly 12 , V. Mainieri 119 , R. Massey 4 , H. J. McCracken 15 - 16 , B. Mobasher 9 , T. 

Murayama 17 , J. A. Peacock 14 , I. Sakelliou 11 , E. Schinnerer 11 , J. D. Silverman 1 , V. Smolcic 11 , Y. Taniguchi 18 , L. 

Tasca 10 , J. E. Taylor 4 , J. R. Trump 8 , G. Zamorani 7 

to appear in COSMOS ApJS special issue 2007 



o 
o 

o 

Q 



> 
o 

m 

(N 
O 

6 ■ 

a: 
. 



ABSTRACT 

We present the results of a search for galaxy clusters in the first 36 XMM-Newton pointings on the 
COSMOS field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 3 x 10~ 15 ergs cm" 2 
s , having one of the widest XMM-Newton contiguous raster surveys, covering an area of 2.1 square 
degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. 
Verification of the cluster candidates is done based on a galaxy concentration analysis in redshift slices 
of thickness of 0.1-0.2 in redshift, using the multi-band photometric catalog of the COSMOS field and 
restricting the search to z < 1.3 and iAB < 25. Wc identify 72 clusters and derive their properties 
based on the X-ray cluster scaling relations. A statistical description of the survey in terms of the 
cumulative log(7V > S) — log(S') distribution compares well with previous results, although yielding 
a somewhat higher number of clusters at similar fluxes. The X-ray luminosity function of COSMOS 
clusters matches well the results of nearby surveys, providing a comparably tight constraint on the 
faint end slope of a = 1.93 ± 0.04. For the probed luminosity range of 8 x 10 42 — 2 x 10 44 ergs s -1 , 
our survey is in agreement with and adds significantly to the existing data on the cluster luminosity 
function at high redshifts and implies no substantial evolution at these luminosities to z = 1.3. 
Subject headings: cosmology: observations — cosmology: large scale structure of universe — cosmol- 
ogy: dark matter — surveys 
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1. INTRODUCTION 

Clusters of galaxies represent a formidable tool for cos- 
mology (e.g. Borgani & Guzzo 2001, Rosati et al. 2002). 
As the largest gravitationally relaxed structures in the 
Universe, their properties arc highly sensitive to the 
physics of cosmic structure formation and to the value 
of fundamental cosmological parameters, specifically the 
normalization of the power spectrum erg and the den- 
sity parameter CIm- Clusters are in principle "simple" 
systems, where the observed properties of the (diffuse) 
baryonic component should be easier to connect to the 
mass of the dark matter halo, compared to the complex- 
ity of the various processes (e.g. star formation and evo- 
lution, stellar and AGN feedback) needed to understand 
the galaxy formation. In particular in the X-ray band, 
where clusters can be defined and recognized as single 
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objects, observable quantities like X-ray luminosity Lx 
and temperature Tx show fairly tight relations with the 
cluster mass (e.g. Evrard et al. 1996, Allen et al. 2001, 
Reiprich & Bohringer 2002, Ettori et al. 2004). Under- 
standing these scaling relations apparently requires more 
ingredients than simple heating by adiabatic compression 
during the growth of fluctuations (Kaiser 1986; Ponman 
et al. 2003, Borgani et al. 2004, 2005). However, their 
very existence and relative tightness provides us with 
a way to measure the mass function (e.g. Bohringer 
et al 2002, Pierpaoli et al. 2003) and power spectrum 
(Schucckcr et al. 2003), via respectively the observed X- 
ray temperature/luminosity functions and the clustering 
of clusters, thus probing directly the cosmological model. 
X-ray based cluster surveys, in addition, can be charac- 
terized by a well-defined selection function, which is an 
important feature when computing cosmological quanti- 
ties as first or second moments of the object distribution. 
Finally, and equally important, clusters at different red- 
shifts provide homogeneous samples of essentially co-eval 
galaxies in a high-density environment, enabling studies 
of the evolution of stellar populations (e.g. Blakeslee et 
al. 2003, Lidman et al. 2004, Mei et al. 2006, Strazzullo 
et al. 2006). 

X-ray surveys in the local ((z) ~ 0.1) Universe, stem- 
ming from the ROSAT All-Sky Survey (RASS, Voges et 
al. 1999) have been able to pinpoint the cluster num- 
ber density to high accuracy. The REFLEX survey, in 
particular, has currently yielded the most accurate mea- 
surement of the X-ray luminosity function (XLF, Ebeling 
et al. 1998; Bohringer et al. 2002; 2004), providing a ro- 
bust z ~ reference frame to which surveys of distant 
clusters can be safely compared in search of evolution. 
Complementary, serendipitous X-ray searches for high- 
rcdshift clusters have been based mostly on the deeper 
pointed images from the ROSAT PSPC (RDCS: Rosati 
et al. 1995; SHARC: Burke et al. 2003; 160 Square De- 
grees: Mullis et al. 2004; WARPS: Perlman et al. 2002) 
and HRI archives (BMW: Moretti et al. 2004); or on 
the high-exposure North Ecliptic Pole area of the RASS 
(NEP: Gioia et al. 2003; Henry et al. 2006). A deeper 
search for massive clusters in RASS (flux limit ~ 10~ 12 
erg cm -2 s -1 ) is carried on by the MACS project (Ebel- 
ing et al. 2001), while the XMM-LSS survey is covering 
~ 9 square degrees to a flux limit ~ 10~ 14 erg cm~ 2 s _1 
(Valtchanov ct al. 2004) 21 . Only recently distant clus- 
ters have started to be identified serendipitously from 
the Chandra (Boschin, 2002) and XMM archives, with 
a record-breaking object recently identified at z = 1.45 
(Stanford et al. 2006). Overall, these results consistently 
show a lack of evolution of the XLF for L < L* ~ 3 • 10 44 
erg s" 1 out to z ~ 0.8, At the same time, however, they 
confirm the early findings from the Einstein Medium Sen- 
sitivity Survey (EMSS, Gioia ct al 1990; Henry et al. 
1992) of a mild evolution of the bright end (Vikhlinin et 
al. 1998a; Nichol et al. 1999, Borgani et al. 2001; Gioia 
et al. 2001; Mullis et al. 2004). In other words, there 
are indications that above z ~ 0.6 one starts finding 
a slow decline in the number of very massive clusters, 

21 All through this paper, we shall adopt a "concordance" cosmo- 
logical model, with H a = 70 km s _1 Mpc -1 , fi M = 0.3, C A = 0.7, 
and — unless specified — quote all X-ray fluxes and luminosities 
in the [0.5-2] keV band and provide the confidence intervals on the 
68% level. 



plausibly indicating that beyond this epoch they were 
still to be assembled from the merging of smaller mass 
units. These results consistently indicate low values for 
f^M and (is ~ 0.7, under very reasonable assumptions on 
the evolution of the Lx — Tx relation. Remarkably, the 
revised value as = 0.76±0.05 yielded by the recent 3-year 
WMAP data (Spergel ct al. 2006), is in close agreement 
to those consistently indicated by all X-ray cluster sur- 
veys over the last five years, both from the evolution of 
the cluster abundance (Borgani ct al. 2001; Pierpaoli 
et al. 2003; Henry 2004) and the combination of local 
abundance and clustering (Schuecker et al. 2003). 

The precision achievable in these calculations essen- 
tially reflects the uncertainty in the relation between 
cluster mass and mcasurablcs such as L x or T x . Progress 
in the knowledge of these relations and their evolution is 
currently limited by the small number of clusters known 
at high redshifts (less than 10 objects known at z > 1), 
although important results have been recently achieved 
(Ettori et al. 2004; Maughan et al. 2006; Kotov & 
Vikhlinin 2005). In addition, virtually all current sta- 
tistical samples of distant clusters have been selected 
serendipitously from sparse archival X-ray images, and 
reach maximum depths of ~ 10~ 14 erg cm~ 2 s . The 
former limits the ability to study distant clusters in the 
context of their surrounding large-scale structure; the 
latter limits to relatively low redshifts the study of low- 
luminosity groups. 

The Cosmic Evolution Survey (COSMOS), covering 2 
□°, is the first HST survey specifically designed to thor- 
oughly probe the evolution of galaxies, AGN and dark 
matter in the context of their cosmic environment (large 
scale structure - LSS). COSMOS provides a good sam- 
pling of LSS, covering all relevant scales out to ~ 50 - 
100 Mpc at z > 0.5 (Scoville et al. 2007a). The rectangle 
bounding all the ACS imaging has lower left and upper 
right corners (RA,DEC J2000) at (150.7988°, 1.5676°) 
and (149.4305°, 2.8937°). To define the LSS, both deep 
photometric multi-band studies of galaxies are carried 
out (Scoville et al. 2007b) as well as an extensive spec- 
troscopy program (Lilly et al. 2007). 

The 1.4 Msec XMM-Newton observations of the COS- 
MOS field (Hasingcr ct al. 2007) provide coverage of 
an area of 2.1 D° to unprecedented depth of 10 -15 ergs 
enr 2 s , previously reached by deep surveys only over 
a much smaller area (Rosati et al. 2002). At the same 
time, the large multi-wavelength coverage makes optical 
identification very efficient (Brusa ct al. 2007; Trump et 
al. 2007). The zCOSMOS program (Lilly ct al. 2007) 
and the targeted follow-up using the Magellan telescope 
(Impcy ct al. 2007) are of particular importance. While 
we anticipate a completion of these programs within the 
next years, we present in this paper the cluster iden- 
tification based on the photometric redshift estimates, 
available as a result of an extensive broad-band imaging 
of the COSMOS field (Capak et al. 2007; Taniguchi et 
al. 2007; Mobasher et al. 2007). 

By design, this paper concentrates on the densest parts 
of the LSS, which have already collapsed to form virial- 
ized dark matter halos of groups and clusters of galaxies, 
populated with evolved galaxies. This approach is quite 
complementary to the study of extended supercluster- 
size structures, which is presented in Scoville et al. 
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(2007b). 

The paper is organized as follows: |JU describes the 
analysis of the XMM-Newton observations of the COS- 
MOS survey; $3] presents our diffuse source detection 
technique; $4] describes the wavelet analysis of the pho- 
tometric galaxy catalog; $5] provides the catalog of the 
identified X-ray groups and clusters of galaxies; iJS] de- 
rives the luminosity function; and |J7] concludes the pa- 
per. 

2. XMM-COSMOS DATA REDUCTION 

For cluster detection, we used the XMM-Newton mo- 
saic image in the 0.5-2 keV band, based on the first 36 
pointings of the XMM-Newton observations of the COS- 
MOS field. The first 25 pointings completely cover the 2 
square degree area of COSMOS at 57% of its final depth 
(Hasinger et al. 2007). 11 pointings have already been 
obtained as a part of the second year observations and 
are included in the present analysis. A description of 
the XMM-Newton observatory is given by Jansen et al. 
(2001). In this paper we use the data collected by the Eu- 
ropean Photon Imaging Cameras (EPIC): the pn-CCD 
camera (Striider et al. 2001) and the MOS-CCD cameras 
(Turner et al. 2001). All Epic-pn observations have been 
performed using the Thin filter, while both Epic-MOS 
cameras used the Medium filter. An analysis of absolute 
frame registration for the XMM COSMOS survey has 
been carried out by Cappelluti et al. (2007). The offsets 
do not exceed 3" and do not affect the current analysis. 

In addition to the standard data processing of the 
EPIC data, which was done using XMMSAS version 6.5 
(Watson et al. 2001; Kirsch et al. 2004; Saxton et al. 
2005), we perform a more conservative removal of time 
intervals affected by solar flares, following the procedure 
described in Zhang et al. (2004). In order to increase 
our capability of detecting extended, low surface bright- 
ness features, we have developed a sophisticated back- 
ground removal technique, which we coin as 'quadruple 
background subtraction', referring to the four following 
steps: 

First, we remove from the image accumulation the pho- 
tons in the energy band corresponding to the Al K a 
line for pn and both MOS detectors and the Si K a line 
for both MOS detectors. The resulting countrate-to-flux 
conversion in the 0.5-2 keV band excluding the lines is 
1.59 x 10~ 12 for pn and 5.41 x 10~ 12 for each MOS de- 
tector, calculated for the source spectrum, corresponding 
to the APEC (Smith et al. 2001) model for a collisional 
plasma of 2 keV temperature, 1/3 solar abundance and 
a redshift of 0.2. 

Second, we subtract the out-of-time events (OOTE) 
for pn. In order to do so, we generate an additional 
event file for each pn event file, which entirely consists of 
events emulating the OOTE, produced using the XMM- 
SAS task epchain. We then apply the same selection of 
events using the good time intervals, generated during 
the light curve cleaning of the main event file. The last 
step is to extract the images with the same selection cri- 
teria (energy, flag, event pattern), normalize them by the 
fraction of OOTE (0.0629191 for the Full Frame readout 
mode used) and to subtract them from the main image. 

Third and fourth, we use two known templates for 
the instrumental (unvignctted) and the sky (vignetted) 
background (Lumb et al. 2002, 2003). The instrumen- 



tal background is caused by the energetic particles and 
has a uniform distribution over the detector. The sky 
background consists of the foreground galactic emission 
as well as unresolved X-ray background and its angular 
distribution on scales of individual observation could be 
considered as flat on the sky, therefore following the sen- 
sitivity map of the instrument (Lumb et al. 2002). To 
calculate the normalization for each template, we first 
perform a wavelet reconstruction (Vikhlinin et al. 1998b) 
of the image without a sophisticated background subtrac- 
tion. We excise the area of the detector, where we detect 
flux on any wavelet scale and then split the residual area 
into two equal detector parts with higher and lower val- 
ues of vignetting than the median value. By imposing 
a criterion that the weighted sum of the two templates 
should reproduce the counts in both areas, we obtain a 
system of two linear equations for two weighting coeffi- 
cients. By solving the latter we derive the best estimate 
for the normalizations of both templates. 

After the background has been estimated for each ob- 
servation and each instrument separately, we produce 
the final mosaic of cleaned images and correct it for the 
mosaic of the exposure maps in which we account for 
differences in sensitivity between pn & MOS detectors 
stated above. The detailed treatment of the background 
is newly developed and has been adopted for all XMM- 
COSMOS projects (Cappelluti et al. 2007). The proce- 
dure of mosaicing of XMM-Newton data has been used 
previously and is described in detail in Briel et al. (2004). 

The resulting signal-to-noise ratio image is shown in 
Fig. [TJ As can be seen from the figure, the image ex- 
hibits a fairly uniform signal-to-noise ratio. Without the 
refined background subtraction, the signal-to-noise im- 
age exhibited large-scale variations, which could mimic 
an extended source. On the image, the circles show the 
position and the angular extent of detected clusters of 
galaxies. The pixel size of the X-ray images employed in 
the current analysis is 4" on a side. 

3. CLUSTER DETECTION TECHNIQUE 

We used the wavelet scale-wise reconstruction of the 
image, described in Vikhlinin et al. (1998b), employ- 
ing angular scales from 8" to 2.1'. We apply the ap- 
proach of Rosati et al. (1998), Vikhlinin et al. (1998b), 
and Moretti et al. (2004) to effectively select clusters 
of galaxies by the spatial extent of their X-ray emission. 
The cluster detection algorithm consists of two parts: (1) 
selection of the area with detectable flux on large angular 
scales, and (2) removal of the area where the flux could 
be explained by contamination from embedded point-like 
sources, most of which are background AGNs, according 
to their optical identification (Brusa et al. 2007). When 
applying this approach to XMM-Newton data, the off- 
axis behavior of the point-spread function (PSF) (e.g. 
Valtchanov et al. 2001) needs to be taken into account. 
The major change in the two-dimensional shape of the 
PSF of the XMM telescope is a tangential elongation 
with respect to the optical axis of the telescope (Lumb 
et al. 2003). In the wavelet scale- wise reconstruction of 
the point source, this produces a flux redistribution be- 
tween the 8" and 16" scales with more flux going into the 
16" at large off-axis angles. At the same time, the flux 
contained in the sum of 8" and 16" wavelet scales does 
not vary much with the off-axis angle and its ratio to 
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Fig. 1. — An image of the signal-to-noise ratio in the 0.5—2 keV band, smoothed with a Gaussian of two pixel width. White color 
corresponds to values smaller than -0.5, while black corresponds to values exceeding 3. Circles indicate the position and size of the detected 
clusters, labeled in correspondence with the cluster number in the catalog. The image is 1.5 degrees on a side. The pixel size is 4" on a 
side. A tangential elongation in the point spread function of XMM-Newton telescopes is clearly visible in the NW corner of the image, 
which has the largest exposure time. 



the flux detected on larger angular scales is greater than 
1.9 at any off-axis angle. Smaller flux ratios have there- 
fore been chosen as a detection criterion for an extended 
source. 

Our wavelet algorithm (Vikhlinin ct al. 1998b) gen- 
erates a noise map against which the flux in the image 
is tested for significance. We assess the sensitivity of 
the survey, by examining the noise map corresponding 
to the 32" scale. The total area over which the source of 
a given flux can be detected is the area where the source 
signal to noise ratio exceeds the chosen detection thresh- 
old of 4. For each cluster this number is unique, yet the 
total flux of the cluster is not recovered at 32". Any 
procedure to recover the total flux introduces a scatter 
in the total flux vs survey area relation. The noise map 
is calculated using the data, and so it includes the ef- 
fect of decreasing sensitivity to a diffuse source detection 
caused by the contamination from other X-ray sources. 
Our calculation of the sensitivity neglects cluster-cluster 



and cluster-point source correlations. In this regard, a 
presence of an unresolved bright cool core could compli- 
cate the detection of a cluster. A modelling of such effect 
will be done within the granted C-COSMOS program. 

The source appearance on the 32" wavelet scale defines 
the area which is further used for measuring the parame- 
ters of the source. The total flux of the source is obtained 
by extrapolation of the flux measured within this area. 
As a first step in selecting the cluster candidates, we clean 
our source list from contamination due to point sources. 
For that, we compare the sum of the two smallest scales 
(8" and 16") with the sum of the largest scales (32", 
64", 128") in the same area to decide whether the area 
contains a diffuse source. If the sum of the two smallest 
scales exceeds the sum of the three largest scales by a 
factor of 1.9 or greater, the source is removed from the 
diffuse source list. This procedure allows us to deal with 
any number of point sources embedded in the same area. 
We find that 80% of the initially selected areas are domi- 
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Fig. 2. — The colors of COSMOS. The wavelet reconstruction of the early-type galaxy concentrations searched in the photo-z catalog 
is color-coded according to the average redshift: blue - 0.2, cyan — 0.4, green - 0.6, yellow — 0.8, red — 1.0. The green contours outline the 
area of the X-ray emission associated with 150 extended source candidates. The image is 1.5 degrees on a side. The pixel size is 10" on a 
side. 



nated by the flux from point sources. Out of the remain- 
ing 150 diffuse source candidates, surviving this test, we 
identify only 76 with peaks in the galaxy distribution, as 
detailed in SjH The bulk of unidentified sources are con- 
fused low-luminosity AGNs, which are not detected on 
the small scales. Some of those sources might be identi- 
fied by chance with a galaxy overdensity. The probabil- 
ity of such event is 1%, as determined by the fraction of 
the total XMM- COSMOS area occupied by the unidenti- 
fied sources, which were initially considered as extended 
source candidates. Sect. 5] describes the construction of 
the catalog of 420 galaxy groups based on the multiband 
photometric data. The number of optical peaks miss- 
ing identification with an extended X-ray source is 345, 



leading to a chance identification for ~ 3 clusters. The 
misidentifications have large separations between the X- 
ray center and the optical center of a group, as a chance 
coincidence of the centers is negligibly small. So, through 
a comparison to optical images we removed 4 sources 
whose shape of X-ray emission was not associated with a 
group of galaxies, finally reducing the cluster sample size 
to 72. For all four rejected systems, one can decompose 
the emission into 3-5 unresolved AGN using the K-band 
images identification, similar to the primary method for 
AGN identification used in the COSMOS field (Brusa ct 
al. 2007). 

It is possible to reduce the number of spurious X-ray 
detections. However, it would result in a loss of survey 
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Fig. 3. — Scale-wise wavelet reconstruction of the XMM-COSMOS mosaic image in the 0.5-2 keV band. Grey color corresponds to the 
surface brightness of 10 -7 counts s" 1 pixel -1 . A surface brightness level exceeding 10 — 5 counts s — 1 pixel -1 is shown in black. Contours 
indicate the location and strength of galaxy structures, identified in the photo-z catalog. The image is 1.5 degrees on a side. The pixel size 
is 4" on a side. 



sensitivity (Vikhlinin et al. 1998b) and is usually done 
to improve the efficiency of the optical follow-up. Since 
such a follow-up already exists in the COSMOS field, 
wc prefer to keep the high sensitivity achieved by our 
detection method. 

Although some additional flux is detected on larger 
scales, the source confusion becomes large on scales ex- 
ceeding 1' at the flux limits typical of our survey. In 
addition, the use of only a single scale for the detection 
has the advantage of allowing a more straight-forward 
modeling. 

Once a diffuse source is detected, the next step is to es- 
timate its total flux. Unlike in the case of point sources, 
this is complicated by the unknown shape of a diffuse 
source. For bright sources, it is possible to carry out a 



surface brightness analysis and estimate the missing flux 
directly. For example, an analysis of the surface bright- 
ness profile for a bright cluster at z = 0.22 in COSMOS 
using the beta model (Smolcic et al. 2007) results in a 
factor of 1.8 higher flux compared to the total flux within 
the aperture defined by the source detection procedure. 

Studies of nearby clusters of galaxies (e.g. Markevitch 
1998) show that most of the cluster flux is encompassed 
within the radius rsoo (a radius encompassing the matter 
density 500 times the critical), so we adopt this radius 
for the total flux measurements, F(< r 50 o). We have 
therefore adopted a method for estimating rsoo, based 
on the iteratively corrected observed flux 

F(< r 500 ) = C p (z,T)F d (1) 
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with the correction, Cp(z,T), defined iteratively 
through the scaling relations as described in §5.11 and 
using the flux and redshift of the system. To estimate 
the observed flux we use the total counts (Fa) in the area 
defined by the detection algorithm (i? t ) and subtract the 
contribution from the embedded point sources based on 
intensity determined by the 4" scale (i?4») of the wavelet. 
The total flux retained in the area is larger than R^i by 
a factor, determined by the size of the area (a) and the 
XMM PSF (Cpsf(o)), with values in the range 1.1-1.5. 
We work in units of MOS1 counts/s on the mosaic. We 
add pn and MOS counts as observed, while putting more 
effective weight to the pn exposure map, so that 

F d = 5.41xlO~ 12 (i?t-Cp S F(a)x J R 4 ") er gss" 1 cm" 2 (2) 

Also, after removal of cool cores, the scatter in the scal- 
ing relations is found to be moderate (Markevitch et al. 
1998; Finogucnov ct al. 2005a). In our survey, a cool 
core is indistinguishable from an AGN and is therefore 
removed from the flux estimates. 

As has been described above, the area for cluster flux 
extraction has been selected based on a particular (32") 
spatial scale. Thus, only a part of the cluster emission is 
used in the detection. This results in an effective loss of 
sensitivity, as not all the cluster emission is used. How- 
ever, the modeling of the detection becomes simple and 
in particular, an influence of the cool cluster cores on the 
cluster detection statistics is reduced. 

4. USE OF PHOTOMETRIC REDSHIFTS TO 
IDENTIFY THE CLUSTERS 

The identification/confirmation of a bound galaxy sys- 
tem would be most robust through optical/NIR spectro- 
scopic redshift measurement, but this is costly for such a 
large survey. Photometric rcdshifts, when measured with 
sufficient accuracy, represent a viable alternative. The 
COSMOS survey, providing both the photo-z (Mobashcr 
et al. 2007) and in the future the spectral-z (Lilly et 
al. 2007) is an ideal field to understand all the pros and 
cons of the use of photo-z to identify the optical coun- 
terparts to X-ray clusters. Diffuse X-ray emission, asso- 
ciated with a galaxy overdensity, is by itself a proof of 
the virialized nature of the object (Ostriker et al. 1995), 
so a combination of diffuse X-ray source and a photo- 
z galaxy overdensity validates the photo-z source as a 
virialized object. Deep multiband photometric observa- 
tions have been carried out in the COSMOS field using 
CFHT, Subaru, CTIO and KPNO facilities (Capak et 
al. 2007). The combination of depth and bands allows 
us to produce an z-band based photo-z catalog with a la 
uncertainty of the redshift estimate of 0.027(1 + z) for 
galaxies with iab < 25 (Mobashcr ct al. 2007) obtained 
without the luminosity prior. While both z and K-band 
data are used in determining the photometric redshifts, 
the input catalog is based on the z-band images. There- 
fore, the redshift range of the cluster search in this paper 
has to be limited to redshifts below 1.3, after which the 
4000A break moves rcdward of the i-band filter. 

To provide an identification of galaxy overdensities in 
the photo-z space, we select from the original photo-z 
catalog only the galaxies, that (1) are classified in their 
SED (spectral energy distribution) as early-types (ellip- 
ticals, lenticulars and bulge-dominated spirals); (2) have 



high quality photo-z redshift determination (95% con- 
fidence interval for the redshift determination is lower 
than Az = 0.4); (3) are not morphologically classified 
as stellar objects. Given the current quality of photo-z, 
we select redshift slices covering the range < z < 1.3 
with both Az = 0.1 and Az = 0.2. To provide a more 
refined redshift estimate for the identified structures, we 
slide the selection window with a 0.05 step in redshift. 
We add each galaxy as one count and apply wavelet fil- 
tering (Vikhlinin ct al. 1998b) to detect enhancements 
in the galaxy number density on scales ranging from 20" 
to 3' on a confidence level of 4ct with respect to the local 
background. The local background itself is determined 
by both the field galaxies located in the same redshift 
slice as well as galaxies attributed to the slice due to 
a catastrophic failure in the photo-z. The fraction of 
galaxies assigned to structures in our analysis varies from 
15% at a z — 0.3 to 2% at z = 1.2. However, the level 
of the background docs not vary strongly with redshift, 
so the detection of the optical counterpart is primarily 
determined by the sensitivity, required to determine the 
photometric redshift for sufficient number of members. 
As will be discussed below, the X-ray detection limits at 
z > 1 in the COSMOS survey correspond to a massive 
cluster, that typically has a sufficient number of bright 
members, enabling the success of the identification pro- 
cedure. 

The angular scales selected for the analysis match the 
extent of X-ray sources found. While the number of 
galaxies within each cluster is not determined by the 
statistics, our ability to see it over the background is 
statistical with the errors determined by the background 
level. Analysis of galaxy overdensities is not sensitive to 
catastrophic failures in the photo-z catalog, as those sim- 
ply reduce the strength of galaxy concentration. Larger 
structures, selected without the prior on the SED type, 
arc reported in Scoville et al. (2007b) and should be 
understood as the LSS environment within which the 
high-density peaks detected via X-ray emission are lo- 
cated. Making a prior on SED type is technically nec- 
essary to increase the contrast in the photometric data. 
From a theoretical point of view, the use of these specific 
SEDs is justified by the large ages of cluster ellipticals 
and the lack of evolution in their morphological fraction 
with redshift (Smith et al. 2005; Postman et al. 2005; 
Wechsler et al. 2005; see also direct estimates from the 
COSMOS field in Capak et al. 2007 and Guzzo et al. 
2007), where the already tested redshift range reaches 
z = 1.1 (Hashimoto ct al. 2005). The ultra-deep field 
of the UKIRT Infrared Deep Sky Survey (Lawrence et 
al. 2006) finds that the galaxy red sequence disappears 
beyond z — 1.5 (Cirasuolo et al. 2006). Moreover, lo- 
cal spiral-rich groups also do not reveal X-ray emission 
associated with their IGM (Mulchaey ct al. 1996). In 
§£] we verify the completeness of the source identifica- 
tion using the V/V max test and conclude that we do not 
lack high-rcdshift identification in our analysis. In com- 
parison with the widely adopted cluster red sequence 
method, the use of photo-z data allows us to include 
fainter members, which is of particular importance for 
galaxy groups. 
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In Fig. [2] we overlay the X-ray contours of diffuse 
sources over the color-coded image of the structures iden- 
tified via the photo-z galaxy selection. The brightness of 
the color is proportional to the number density of galax- 
ies, while the color represents the average redshift. In ad- 
dition, to show the correspondence between galaxy struc- 
tures and X-ray emission, we provide in Fig. [3] an overlay 
of the wavelet-reconstructed X-ray image with the con- 
tours, showing the sum of the wavelet-filtered redshift 
galaxy slices. 

To determine the redshift of an X-ray structure, we cal- 
culate the position and maximum of the galaxy density 
peak using the photo-z slices. The extended source can- 
didate is considered identified as a cluster if it contains a 
galaxy density peak inside the detected extend, which is 
typically 20" — 40". The redshift and center of the clus- 
ter is selected to be associated with the strongest galaxy 
peak inside the cluster. Due to the wings in the redshift 
estimate using the photometric technique, the cluster is 
detected over a range of redshifts. Selecting the peak 
value is thought to yield the most likely redshift of the 
cluster. In the case of several overlapping structures, this 
could lead to biases due to redshift selection effects and 
in any case should be verified. In one instance (cluster 
133), we manually put a preference for a high redshift 
counterpart, which matched both the X-ray center and 
the redshift of BCG of the system. In total, 72 X-ray 
clusters have been identified using this method. Their 
properties are listed in Tab. 1 and discussed below. 

5. A CATALOG OF IDENTIFIED X-RAY 
CLUSTERS 

In this section we describe our catalog of 72 X-ray 
galaxy clusters detected so far in the COSMOS field. In 
the catalog we provide the cluster identification number 
(column 1), R.A. and Deck of the peak of the galaxy 
concentration identified with the extended X-ray source 
in Equinox J2000.0 (2-3), the estimated radius r$oo in 
arcminutes (4), the cluster flux in the 0.5-2 keV band 
within rsoo in units of 10~ 14 ergs cm -2 s _1 with the cor- 
responding 1 sigma errors (5), photometric redshift (6), 
an estimate of the cluster mass, M500, in units of 10 13 M Q 
(7), rest-frame luminosity in the 0.1-2.4 keV band (8), an 
estimate of the IGM temperature in keV (9), the source 
of redshift (10) with 1 denoting the photometric redshift 
estimates and denoting the redshift refinement using 
information available from the SDSS DR4 (York et al. 
2000, Abazajian et al. 2005; Adelman-McCarthy et al. 
2006), the results of the pilot zCOSMOS program, and 
the first results of the main zCOSMOS program (Lilly et 
al. 2007). The refinement of the photo-z is complete for 
z < 0.2, which reduces the degree of uncertainty in calcu- 
lation of distance-dependent properties of the sample to 
a level, well within the scatter in the scaling relations for 
X-ray clusters (e.g. Stanek et al. 2005). The formal error 
on the positional uncertainty is of the order of 10". The 
errors provided on the derived properties are only statis- 
tical. A future study of systematic errors will be done by 
direct comparison between our reported mass estimates 
and the weak lensing measurements (Leauthaud et al., 
in prep.; Taylor et al., in prep.). For a number of clus- 
ters with best X-ray statistics, it appears that the peak 
of the galaxy density can be somewhat offset from the 
X-ray center, which in turn has a better match with the 



position of the brightest cluster galaxy. Further investi- 
gation of this interesting property will be done elsewhere. 

In order to calculate the properties of the clusters, we 
use the known scaling relations at low-rcdshift and evolve 
them according to the recent Chandra and XMM ob- 
servations of high-rcdshift clusters of galaxies (Kotov & 
Vikhlinin 2005; Finoguenov et al. 2005b; Ettori et al. 
2004; Maughan et al. 2006). Following Evrard et al. 
(2002), we choose not to correct for the redshift evolu- 
tion of the ratio of the density encompassed by the cluster 
virial radius to the critical density. 

5.1. Flux estimates 

Even a basic statistical description of the cluster sam- 
ple, such as log(A) — log(S'), is not entirely straight- 
forward, as it requires a knowledge of the level of the 
emission at large distances from the cluster center. Most 
cluster surveys to date either use much deeper reobscr- 
vation of the clusters for this purpose or extrapolate the 
detected cluster flux assuming some model for the spatial 
shape of the emission. Although most observers agree 
that the /3-model (Jones & Forman 1984; 1999) is a good 
description of the cluster shape, characterization of the 
outskirts of groups and clusters of galaxies is still uncer- 
tain (e.g. Vikhlinin et al. 2006; Borgani et al. 2004). 
Moreover, groups reveal a large scatter in the shape of 
their emission (Finoguenov et al. 2001; Mahdavi et al. 
2005). 

For the purpose of this paper, we extrapolate the mea- 
sured flux (Fd) by an amount of Cp(z,T), which uses 
the /3-model characterization of cluster emission, with 
the following parameters 

= 0.4 x (fcT/lkeV) 1/3 (3) 

and 

rcore = 0.07 x (fcT/lkeV) 63 x r 50 o (4) 

We note that although the cool core clusters require 
a second spatial component to describe the central peak 
of their emission, in our case the contribution of a cool 
core cannot be distinguished from point sources and has 
therefore been removed from the flux estimates. Since 
only the large-scale component exhibits a scaling with 
temperature (Markevitch 1998), removal of the center 
should result in reduced scatter in the derived cluster 
characteristics, such as the total mass. More impor- 
tantly, exclusion of cool cores from the detection pro- 
cedure decreases an observational bias towards low-mass 
groups with strong cool cores, and make a selection of 
the sample to be closer to mass selection, and reduces the 
bias discussed in Stanek et al. (2006) and O'Hara et al. 
(2006), while facilitating the modelling of the cluster de- 
tection. However, subtle differences in the derived char- 
acteristics of the sample are expected as discussed below. 
In calculating the rest-frame luminosity, we iteratively 
use the total flux within an estimated r$oo (Cp(z, T)Fd) 
and apply the K-correction (K(z, T)) accounting for the 
temperature and redshift of the source, following the ap- 
proach described in Bohringer et al. (2004) and assum- 
ing an element abundance of 1/3 solar. So, finally, the 
derived luminosity is 

£o.i-2.4kcV = ^d 2 L K{z, T)C p (z, T)F d (5) 



10 



Finoguenov et al. 




13— « 



10-' 5 10-' 4 10-' 3 0.1 0.5 1 5 10 
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Fig. 4. — Left panel: Cumulative cluster number counts (log(iV > S) — log(S')) for the COSMOS field. The solid histogram shows the 
data and grey histograms denote the 68% confidence interval. The black solid/short-dashed curve shows the results of the modeling of 
RDCS (Rosati et al. 2002), with the solid part corresponding to fluxes sampled by their data, while the short dashed line denotes the 
model prediction. The long dashed curve shows the prediction for no evolution in the luminosity function in Rosati et al. (2002), which 
provides a good fit to our data. The dashed histogram shows a typical difference due to assumption of the scaling relations. Right panel: 
Survey area as a function of the total source flux in the 0.5-2 keV band. Open circles show the area corresponding to the source, while the 
solid line shows the sensitivity curve for the flux in the detection cell, scaled up to match the displayed data. The scatter of points around 
the curve is determined by differences in the flux losses between distant and nearby sources. 



Wc note that a similar approach for reconstructing the 
cluster flux is taken in Henry et al. (1992), Nichol et 
al. (1997), Scharf et al. (1997), Ebeling et al. (1998), 
Vikhlinin et al. (1998b), de Grandi et al. (1999), Reprich 
& Bohringer (2002), and Bohringer et al. (2004). A 
choice of r 500 as a limiting radius is motivated by the 
results of Vikhlinin et al. (2006) on the steepening in 
the surface brightness profiles, observed near this radius. 

To estimate the temperature of each cluster, we use 
the I/o.i-2.4keV — T relations of Markevitch (1998) for 
the case of excised cool cores: 

kT = 6 keV X l ( 1 °gio(^o.i-2.4 k c V £r 1 )-44.45)/2.1) (g) 

where 

e z = (n M (i + zf + n A y/ 2 (7) 

The estimates of the total gravitational mass and the 
corresponding rsoo are performed using the M-T rela- 
tion (F. Pacaud 2005, private communication) re-derived 
from Finoguenov et al. (2001) using an orthogonal re- 
gression and correcting the masses to hro and a ACDM 
cosmology: 

M 500 = 2.36 x 1O 13 M x T 189 ^" 1 (8) 

and 

r 500 = 0.391 Mpc x (fcT/keV) 63 ^" 1 (9) 

Although the L 500 - T relation of Markevitch (1998) is 
formally derived for high-temperature clusters, a com- 
parison of that relation with the new results on the 
^500 ~ T relation for groups (Ponman T.J., private com- 
munication), shows that the two are very similar. Based 
on direct comparison between the predicted and mea- 
sured temperatures for a number of clusters in the COS- 
MOS (e.g. Smolcic et al. 2007; Guzzo et al. 2007), the 
mass estimates provided in Tab. 1 should be good to a 
factor of 1.4, due to both the scatter in the scaling rela- 
tions and uncertainty in our knowledge of their redshift 
evolution. 



5.2. Cluster counts 

It is common to characterize a cluster survey by its 
area as a function of the limiting flux and present the 
results as a relation between a cumulative surface density 
of clusters above a given flux limit vs the flux value, the 
cluster log(A > S) - log(S) (e.g. Rosati et al. 1998). 
In order to take into account the difference between the 
total flux and the observed flux, the log (AT > S) — log(<S) 
is computed by assigning to each cluster a weight equal 
to the inverse value of the area for its observed flux. Such 
area-flux relation is determined by the wavelet detection 
algorithm and is therefore known precisely. This weight 
is further added to the flux bin in correspondence to the 
total flux of the cluster. 

While the calculation of the survey area as a function 
of cluster flux in the detection cell is exact, the survey 
area as a function of total cluster flux is known only 
approximately and exhibits a scatter, as illustrated in 
the right panel of Fig. [4] Lower redshift systems require 
more extrapolation for a given flux in the detection cell. 
Similarly, within the same redshift range, more massive 
systems require larger degree of flux extrapolation given 
the fixed size of the detection cell. 

The left panel of Fig. H shows the log(A > S) - log(5) 
relation of the COSMOS field clusters. Although with a 
somewhat higher normalization, the COSMOS log(iV > 
S) — log(S') is statistically consistent with RDCS results 
of Rosati et al. (2002) for the fluxes S > 10" 14 ergs 
cm~ 2 s" 1 . While a similarly higher normalization of the 
log N — log S has also been reported for the XMM-LSS 
survey (Pierre et al. 2005) as well as for the 160 square 
degrees survey (Vikhlinin et al. 1998b), an important 
difference between those surveys and the RDCS consists 
in extrapolation of the cluster flux beyond the detection 
radius. As no extrapolation has been done for RDCS, 
the difference in the results is due to higher flux being 
assigned to each source, and not due to the higher source 
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density. At fluxes below 10 -14 ergs cm" 2 s" 1 the XMM- 
COSMOS is the first survey to yield rich observational 
data, allowing us to determine the log(A > S) — log(<S) 
with good statistics down to S ~ 10" 15 ergs cm" 2 s _1 . 
We note that the prediction for no evolution in the lumi- 
nosity function obtained by local surveys, as summarized 
in Rosati et al. (2002), provides a good fit to our cluster 
counts. 

To check the uncertainties in our derived \og(N > 
S) — log(S') relation due to our procedure for the esti- 
mate of the total flux, in Fig. [4] we show also the counts 
obtained if we compute the total flux for each clusters 
assuming /3 = 0.6, r corc = 0.08r5oo and using the param- 
eters of the M — T relation from Vikhlinin et al. (2006). 
The number counts derived with these assumptions are 
close to the upper envelope of the 1 sigma confidence re- 
gion of our \og(N > S) — logos'). This difference in nor- 
malization is due to the fact that the M — T of Vikhlinin 
implies a larger size of the clusters, in particular at the 
lower mass end, compared to the relation adopted here. 
We would like to point out several caveats in the use 
of the presented \og{N > S) — log(S'). One is that clus- 
ter identification at z > 1.3 has not been included, which 
can increase slightly the faint end counts. Another caveat 
consists in the application of the log( A > S) — log(S) for 
predicting the number of clusters in surveys, as it should 
take into account both the relation between the total flux 
and the flux of the cluster in the detection cell, which is 
sensitive to the adopted detection method. Finally, re- 
moval of the cool cores from the flux estimates, results 
in an underestimate of the total X-ray flux from clusters, 
typically by 20% based on cluster studies at intermediate 
rcdshifts (e.g. Zhang et al. 2006). A study of influence of 
cool cores on both log(A > S) — log(S) and the evolution 
of X-ray luminosity function will be carried out within 
the granted Chandra program. 

5.3. Sample characteristics 

In Fig. [3] we plot the observed characteristics of the 
XMM- COSMOS cluster sample together with detection 
limits implied by both survey depth and our approach 
to search for clusters of galaxies. The solid gray line has 
been calculated by requiring that the 2 x r 50 o = 32" and 
shows which sources cannot be detected as extended by 
our technique. The dotted grey line is calculated im- 
posing a criterion 2 x r core = 32" and shows for which 
clusters the cores could be resolved in our method. The 
dashed grey line is calculated by requiring that r.500 = 
32" to reveal clusters for which removal of the embed- 
ded point sources will result in a strong underestimate 
of the cluster emission. The black lines show the detec- 
tion limits of the survey achieved over 90, 50 and 10% of 
the total area. This comparison shows that with a given 
method it is possible to go a factor of 10 deeper without 
losing X-ray groups of galaxies that would appear point- 
like, which explains the success of the application of our 
cluster detection method to deep fields. Thus, a use of 
32" scale does not introduce a systematics towards the 
source detection, as for the flux limits of the survey the 
expected extent of the X-ray emission is larger. On the 
other hand, resolving the cores for the detected groups 
and clusters is difficult and requires a more refined PSF 
modeling, that we postpone to a future effort. In all 
cases, the central 0.5rsoo is resolved. An underestimate 



of mass due to oversubtraction of the unresolved cluster 
emission would occur between the long-dashed and solid 
grey lines, and is currently below the survey sensitivity 
limits. We have already demonstrated that clusters with 
rsoo smaller than our detection cell cannot be accessed 
based on the sensitivity limits and do not introduce any 
systematics in the selection. Likewise, the clusters much 
larger than the detection cell, which are the brightest at 
each redshift, are all detected de facto, because the avail- 
able statistics allows us to map the core of the cluster at 
32" resolution. In a much shallower survey, where the 
detection limits would be similar to the limiting flux at 
which the cluster cores are resolved (e.g. for XMM it is 
the short-dashed curve in Fig. 5), a problem with a fixed 
detection size might occur. Thus we conclude, that the 
success of the method is a result of a good match be- 
tween the detection cell, the depth of the survey and the 
properties of X-ray emission of clusters of galaxies. 
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Fig. 5. — Illustration of the cluster luminosity probed as a func- 
tion of cluster redshift in the XMM COSMOS survey. Filled circles 
represent the detected clusters with error bars based on the sta- 
tistical errors in the flux measurements only. Short-dashed, long- 
dashed and solid black lines show the flux detection limits associ- 
ated with 90, 50 and 10% of the total area, respectively. Grey lines 
indicate the limits imposed by the detection method (i.e. due to 
the angular resolution of the XMM-Newton telescopes). Systems 
below the short-dashed grey line have unresolved cores, systems 
below the long-dashed grey line (none in our sample) suffer from 
oversubtraction of core emission. No system below the solid grey 
line can be detected as an extended source using the method de- 
scribed in this paper. 

In Fig. [5] we report the redshift distribution of both 
galaxy overdensities and the identified X-ray structures. 
Galaxy overdensities are determined from the analysis 
described in SjH In particular, wc display the histograms 
for the wavelet detected peaks in the photo-z slices with 
Az = 0.2, which is also used to normalize the abundance 
of galaxy structures in Fig. [6] 

Comparison between all galaxy overdensities with the 
X-ray selected systems shows that 25% of the galaxy sys- 
tems are identified in X-rays at z < 0.5. At higher red- 
shifts this proportion drops to 15%. The reason for a 
change in the identification rate consists in the increase 
of the total mass of a group of galaxies that could be de- 
tected through its X-ray emission in the XMM-COSMOS 
survey. At the same time, the depth of the catalog al- 
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lows us to cover the L* galaxies even at z ~ 1.2 and 
the abundance of galaxy counterparts is not limiting the 
identification process at z < 1.3. 
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Fig. 6. — Differential redshift distribution (dN/dz) of 420 photo- 
z galaxy concentrations (grey histogram) as well as 72 X-ray (black 
histogram) clusters of galaxies in the COSMOS field. 

The depth of the X-ray survey could also be char- 
acterized by the median redshift of the catalog, which 
is shown in Fig. [71 The median redshift for the to- 
tal sample is 0.45. We define a subsample of the 12 
most massive objects in the sample: the galaxy clus- 
ters, which are traditionally defined as M200 > 10 14 M Q , 
where M200 — 1.7Msoo- This subsample is character- 
ized by a median redshift of 0.75 and should not be af- 
fected by our limiting X-ray flux, yet this value could 
be somewhat underestimated, due to our z < 1.3 limita- 
tion. Both the surface density of clusters and the median 
redshift matches well the expectations for the future SZ 
surveys. Thus, COSMOS field can be used to calibrate 
the efficiency of SZ surveys. 
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Fig. 7. — Normalized cumulative (N(> z)/N total ) redshift dis- 
tribution of the diffuse X-ray sources, identified with a galaxy con- 
centration based on the photo-z's. The thick grey line represents 
the full sample of 72 clusters, and the solid black line represents 
the 12 clusters (M 20 o > 1O 14 M ). 



6. X-RAY LUMINOSITY FUNCTION 

The volume probed by the COSMOS survey is rep- 
resentative for LSS studies (Scoville et al. 2007a) and 
the size of the COSMOS X-ray cluster catalog is simi- 
lar to many cluster surveys, which is due to a combina- 
tion of high sensitivity achieved over a sufficiently large 
area. For comparison, the volume probed by the XMM- 
COSMOS survey to a redshift of 1.3 is achieved by RE- 
FLEX (Boehringer et al. 2001) by z ~ 0.05, the 400D° 
survey (Burenin et al. in prep.) by z ~ 0.2, and RDCS 
(Rosati et al. 2002) by z ~ 0.4. At these redshifts, all 
those surveys provide volume-limited samples of clusters 
brighter than a few times 10 43 ergs s _1 . For more lu- 
minous systems, wider surveys provide more volume, so 
the strength of XMM-COSMOS consists in probing the 
redshift evolution of 10 43 ergs s _1 clusters. To provide a 
more detailed comparison we construct here a luminos- 
ity function for the clusters in our sample. The new flux 
regime probed by the COSMOS survey yields an abun- 
dance of low-mass groups, so our results will contribute 
to the refinement in the faint end of the luminosity func- 
tion. 

In defining the luminosity function we followed the ap- 
proach described in Mullis et al. (2004). In particular, we 
used their Eq.4, which is adopted from Page & Carrera 
(2000) on the refined estimate of the luminosity function. 
In calculating the maximum volume probed by the sur- 
vey as a function of cluster luminosity, we first tabulate 
the limiting luminosities on the grid defined by tabula- 
tion of the survey area vs observed flux and sample the 
redshifts 0-1.3 with a step of 0.01. We take into account 
the typical extent to which the source is detected, which 
is found to be a function of the observed flux, and es- 
timate the degree of extrapolation performed (Cp(z, T)) 
in obtaining the total flux and apply the K-correction 
(K(z,T)). 

In the no-evolution case and in the absence of strong 
clustering, one expects to detect the clusters uni- 
formly throughout the probed volume ((V/V max ) = 0.5, 
Schmidt 1969). So, if any issues of incomplete identifi- 
cation, e.g. resulting from the use of the photo-z cat- 
alog, are important at some redshift (e.g. at very low 
redshifts or at very high redshifts), this would introduce 
a distortion in the distribution of clusters over the vol- 
ume. To check if there are biases for some class of objects 
(e.g. low-luminosity objects), we plot in Fig. [5] the ratio 
between the volume towards the system and the max- 
imum volume at which it could be detected. As some 
values exceed unity due to the scatter in the flux-area 
relation, we replaced them by unity for illustrative pur- 
poses. The numbers of clusters above and below 0.5 are 
roughly equal at all luminosities, indicating no large se- 
lection effects. The mean value of V/V max for the survey 
is equal to 0.48 ±0.06, which is consistent with 0.5 within 
the statistical errors. This is an important result in itself, 
that may illustrate that morphological changes observed 
in high-redshift clusters (Postman et al. 2005) do not 
cause strong rcdshift-dcpendcnt selection effects. It is 
clear, however, that any effects of incompleteness occur- 
ring at the 10% level would be hard to detect with the 
size of our sample. 
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Fig. 8. — Test for the sample redshift completeness (V/V ma x). 
The estimates exceeding 1 (due to the scatter in the flux-area re- 
lation) are substituted with 1. The number of clusters above and 
below 0.5 is roughly equal at all luminosities, indicating no large 
selection effects. 

Finally, in Fig. [9] we present the luminosity function of 
XMM COSMOS clusters. A comparison with the results 
of the REFLEX (Boehringer et al. 2001) and the BCS 
survey (Ebeling et al. 1997) displayed in Fig. [9] shows 
that, at the luminosity range probed by the COSMOS 
survey, the evolution in the luminosity function is not 
statistically significant. 

The COSMOS data allow us to put tight constraints on 
the slope of the faint end of the luminosity function. To 
characterize it, we fitted a Schechter function to the data 
adopting the L* x and 4>* parameters in correspondence to 
the best fit values of BCS survey (9.1 x 10 44 /i^ 2 ergs/s _1 
and 7.74 x 10 _8 /if Mpc~ 3 ). We achieve an acceptable 
value of the reduced \ 2 — 0.7 for 7 degrees of freedom and 
constrain the value of the slope to a = 1.93±0.04 (where 
dn/dL oc L~ a ), using the Gchrcls (1986) approximations 
in calculating the confidence limits for the case of small 
number statistics. The value of the slope compares well 
with the BCS result of a — 1.85 ± 0.09. Our slope value 
is also within the uncertainty reported for the 160 square 
degrees survey (Mullis et al. 2004). However, their use of 
the 0.5 — 2 kcV energy band results in a somewhat lower 
value of the slope, and a strict comparison of luminosity 
functions is difficult. The lower number of groups in 
the REFLEX survey (dotted line in Fig. [9]) is thought 
to be due to a combination of the small survey depth 
at low luminosities and a presence of a local southern 
void, where most of the survey area is located and whose 
effect has been demonstrated through differences within 
the sample (Bohringer et al. 2002). 

To illustrate the lack of redshift evolution in the lu- 
minosity function, in Fig. [9] we split the sample in two 
redshift bins, 0-0.6 (dotted crosses) and 0.6-1.3 (solid 
crosses), retaining only the luminosity bins derived us- 
ing at least three clusters. The two sub-samples overlap 
only in a single luminosity bin, where the corresponding 
cluster abundances agree within the errors. Since the 
high- luminosities are well probed only at rcdshifts higher 
than 0.6, the good match between our measurements and 



the local luminosity function is an indication of the ab- 
sence of a significant redshift evolution. This finding is 
in agreement with the results of Mullis et al. (2004), 
where detectable evolutionary effects are seen just above 
L x ~ 10 44 ergs/s. We note that our measurements are 
of comparable quality to the existing data compiled in 
Mullis et al. (2004) and provide a refinement to the 
knowledge on the luminosity function at high redshifts. 
A further refinement of the COSMOS survey stems from 
its ability to measure the mass independently via weak 
lcnsing and therefore reduce the effects of systematics in 
determining the mass function. The study of cluster evo- 
lution provides competitive constraints in the — 
plane, and yields results in agreement with other mea- 
surements (Vikhlinin et al. 2003). Finally, we note that 
both our cluster counts and the luminosity function are 
consistent with no evolution in the luminosity function 
in the 8 x 10 42 — 2 x 10 44 ergs s _1 range. This provides 
further evidence in favor of a consistent modeling of the 
XMM-COSMOS survey sensitivity presented in this pa- 
per. 
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Fig. 9. — Luminosity function of clusters in the COSMOS field. 
Dotted crosses indicate the data in the redshift range 0-0.6, grey 
points are the data in the redshift range 0-1.3 and solid crosses 
indicate the data in the redshift range 0.6—1.3. The dotted line 
shows the luminosity function of the REFLEX survey (0 < z < 
0.3, Bohringer et al.2001) and the dashed line shows the results 
of BCS survey (Ebeling et al. 1997), which illustrates the current 
uncertainly on the shape of the luminosity function at z < 0.3. 
The solid line shows the best fit to the COSMOS data. 



7. SUMMARY 

We present a description of our X-ray based cluster de- 
tection method and the first results of the cluster search 
using the XMM-COSMOS survey. Our flux range is 
3 x 10" 15 - 10~ 13 erg cm" 2 s" 1 in the 0.5-2 kcV band. 
We run a separate analysis of the photo-z catalog to iden- 
tify 420 early- type galaxy concentrations, which provide 
an identification to 72 X-ray cluster candidates. We fur- 
ther present the statistics for those clusters in terms of 
log A — log 5, dN/dz and dn/dL. By comparison with 
local cluster surveys, we find no evolution in cluster num- 
ber abundance out to a redshift of 1.3 in the luminosity 
range of L . i~2. 4kcV : 8x 10 42 — 2 x 10 44 ergs s _1 . This fur- 
ther implies that the surface density of clusters detected 
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in the flux range 10~ 15 — 10~ 14 erg cm" 2 s" 1 should cor- 
respond to the prediction of no evolution, higher than 
implied by Rosati et al. (2002). Such high surface den- 
sity of clusters has been found by both COSMOS and 
XMM-LSS surveys (Pierre et al. 2005). The published 
results on deeper X-ray surveys, CDFS and CDFN, con- 
tradict this view (Rosati et al. 2002). It is therefore 
important to understand the origin of this inconsistency, 
which is likely related to evolution of cluster X-ray lumi- 
nosity function at z > 1.2, probed at fluxes fainter than 
10~ 15 erg cm~ 2 s _1 . 

X-ray selected clusters provide well-defined informa- 
tion on the most massive high density peaks in the COS- 
MOS field and allow a follow-up study of the evolution 
in the galaxy morphology in dense environments. 
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